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CHAPTER  I 

GENERAL  DEVELOPMENTS 

This  thesis  is  a  study  o£  the  demodulation  of  antipodal  signals 
transmitted  through  a  nonlinearly  distorting  channel.  The  motivation  for 
this  work  is  the  related  research  effort  by  Poor,  reported  in  [1]  and  [2], 
General  expressions  will  be  drawn  freely  from  this  literature. 

An  antipodal  signalling  problem  involves  the  discrimination  between  two 
mutually  exclusive  events.  To  illustrate,  consider  the  following  pair  of 
hypotheses  for  an  observation  Z  a  (z(t);  0  ^  tS  t}: 

HQ:Z(t)  -  Tl(t)  -  s(t)  ;  0  £  t  <  T 

versus 

H1:Z(t)  -  7l(t)  +  s(t)  ;  0<tSI,  (1.1) 

where  (s(t);  0  <  t  <  t}  is  the  known  deterministic  signal,  and 
[T|(t);  0  s  t  <  l]  is  a  sample  function  of  a  zero-mean  additive  Gaussian 
noise  process. 

To  discriminate  between  Hq  and  in  (1.1),  consider  the  decision  rule: 

'  1  ;  if  <h,z>  a  0 

y(h;z)  =• 

.-1  ;  if  <h,z>  <  0  ,  (1.2) 

where  h  s  £h(t);  0<  t  <  T)  denotes  the  impulse  response  of  a  linear  time- 
invariant  filter,  and  the  operation  (h,z)  is  given  by 

T 

<h,z>  -  f  h(T-t)z(t)dt  . 

J0 


(1.3) 


2 


For  a  value  of  Y(h;z)  ■  1,  the  decision  rule  indicates  the  choice  H^, 
while  Y(^;z)  *  -1  will  indicate  the  choice  Hq.  Assuming  Hq  and  are 
equally  likely  to  occur,  the  probability  of  error  using  decision  rule 
Y(h;z)  is  given  by  the  equation 

pe(h)  »  l-5C<h,s>/[<h,RNh>]^},  (1.4) 

where  §{x)  is  the  unit  normal  cumulative  distribution  function 

t  X  , 

$fx)  -  (2tt)  *  j  exp[-u4/2]du  ,  (1.5) 

-0» 

(h,s)  is  as  in  (1.3),  and 

T  T 

<h,R  h>  *  J*  h(T-t)  J*  RYt,u)h(T-u)du  dt,  (1.6) 

0  0 

where  s  {^(t.u)  ;  0  <  t;  u  <  t}  is  the  known  autocorrelation  function 
of  the  additive  noise  process  Cn(t);  0  ^  t  <  t}  . 

The  probability  of  error  pg(h)  Is  a  primary  measure  of  system 
performance  for  digital  data  communications .  For  optimal  design  with 
respect  to  this  measure,  signal  detection  is  assumed  to  be  accomplished 
by  a  linear  receiver.  Synchronous,  or  coherent,  detection  in  an  additive 
Gaussian  noise  background  requires  a  matched-filter,  or  equivalent 
correlation  detection,  to  achieve  minimum  probability  of  error  for  fixed 
signal  and  noise  conditions.  The  matched-filter  impulse  response  h(t) 
is  known  to  be  the  solution  to  the  integral  equation 

T 

s(t)«J*  RTJ(t,T)h(T-T)dT  ,  (1.7) 

0  " 

a  Fredholm  equation  of  the  first  kind.  Here  the  functions  s(t)  and 
R^(t,r)  are  as  previously  described. 


The  choice  of  antipodal  signalling  also  minimizes  probability  of  error. 


since  this  signal  design  becomes  optimal  for  general  binary  receivers  by 
maximizing  signal- to-noise  ratio  (SNR).  The  signal  Interval  and  the 
receiver  observation  interval  are  also  assumed  to  coincide. 

The  nonlinear  channel  distortion  is  modelled  to  affect  the  signal 
structure  alone,  not  the  noise  covariance.  Examples  of  this  phenomenon 
include  transmitter  or  receiver  generated  noise,  improper  phase  locking, 
channel  fading,  and  multiple  transmission  paths.  It  is  usually  convenient 
to  model  all  distortions  in  the  communications  channel,  and  to  assume  an 
ideal  transmitter  and  receiver.  This  channel  distortion  motivates  the 
following  model  for  the  received  signal  class  J  as  found  in  the  work  by 
Poor  [ 1 ] ; 


J  31  [s  €  L2[0,T]:  JMs(t)  -  s° ( t )  | 2  dt<  A}, 


U.8) 


where  {s^(t);  0  ^  t  <  T]  is  a  known  nominal  signal  model,  L2[0,T]  denotes 
the  square-integrable  functions  on  the  interval  [0,T],  and  A  is  the  degree 
of  channel  distortion. 

To  design  a  suitable  filter  for  detection  within  the  distortion  model  of 
(1.8),  we  consider  a  maximin  design  criteria  as  in  [1],  i.e.,  we  seek  a 
solution  to 

max  min  t<h,s>/[<h,ilh>]%}  .  (1.9) 

h  s  €  J  M 

A  solution  to  (1.9)  will  have  the  largest  possible  minimum  output  signal- to- 
nolse  ratio  and  thus  its  worst-case  performance  will  be  the  best  possible 
within  the  tolerances  set  by  the  uncertainties  in  the  signal  structure. 

Such  a  filter  is  known  as  a  robust,  or  optimally  stable,  matched  filter. 


«* 


Given  Che  problem  of  interest,  with  a  known  noise  autocovariance 

function  [r^}  and  a  signal  structure  as  given  in  (1.8),  the  robust  filter 

£ 

solution  is  that  impulse  response  h  (t)  which  is  the  solution  to  the 
equation 

s°  *  (R^  +  aQ  I)hR  .  (1.10) 

this  is  a  general  result  as  given  in  [2],  where  s^  is  the  nominal  signal 
waveform,  1^  denotes  the  identity  operator  on  L^fO.T],  and  cr^  is  a  non¬ 
negative  constant  determined  by  the  equation 

T 

a2Q  J*JhR(t)|2  dt  =•  A  .  (1.11) 

Further,  it  is  shown  that  the  impulse  response  h  (t)  is  the  matched- filter 

2 

corresponding  (or  matched)  to  the  least  favorable  signal  s  (t)  €  J  given  by 

sL(t)  *  s°(t)  -CTQhR(t)  ,  (1.12) 

and  it  follows  in  [2]  that  the  worst  filter  performance  over  the  class  of 
signals  J  and  noise  is 

max  p  (hR(t)) *  1  -  ${[  nin  <hR,s)/[<hR,R^hR> ]^J 
s(t)€  J  6  s(t)€*' 

*  l-$C[<sL,hR>]*]  .  (1.13) 

Note  that  (1.10)  specifies  the  robust  filter,  and  (1.13)  specifies  the  worst 
performance  of  this  filter.  Poor  £1-1  points  out  that  since  the  identity 
operator  £  corresponds  to  white  noise,  the  effect  of  distortion  of  the  type 
of  (1.10)  is  equivalent  to  the  effect  of  adding  white  noise  of  spectral 
height  Oq  to  the  communications  channel. 
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Equation  (1.10)  is  written  explicitly  as 


T 

8°(t)  -  aQhR(T-t)  +  J*  Rj,( 


(1.14) 


which  is  a  Fredholm  equation  of  the  second  kind.  Properties  of  this 

equation  as  well  as  solution  techniques  are  found  in  the  literature  [3]. 

This  thesis  will  employ  the  method  of  Hilbert  and  Schmidt.  In  particular, 

2 

if  R^(t,u)  is  continuous  on  [0,T]  ,  then  a  continuous,  square-lntegrable 

D 

solution  for  h  (t)  will  always  exist. 

2 

Moreover,  R^t.u)  has  a  Mercer  expansion  on  [0,T]  given  by 

00 

^(t.u)  -  Z  \,tN(t)tN(u)  ,  (1.15) 

2 

with  uniform  convergence  on  [0,T]  .  The  constants  f\„;  N  -  1,2,3,...} 

N 

and  the  functions  0{>N(t);  n  *  1,2,3,...}  are  the  eigenvalues  and  eigenfunctions, 
respectively,  of  [l^} ;  they  are  the  solutions  to  the  homogeneous  equation 

T 

X*(t)  -  J  RJl(t,u)t(u)du;  0  <  t  <  T  ,  (1.16) 

with  the  tN(<0's  being  orthonormal;  i.e.. 


I0  ♦h(u)*m<u>‘1u 


(  1 


\  o 


N  -  M 


N  t  M  . 


(1.17) 


It  is  desirable  to  obtain  series  solutions  to  the  quantities  of 

Interest  which  will  represent  a  unique  solution  to  (1.14).  In  Lovitt  [3], 

£ 

it  is  shown  that  the  unique  solution  h  (t)  is  given  by  the  series 

hR(T-t)  -  a0[s°(t)  -  2  cN(l  +  o0/XN)"SN(t)];  0  <  t  <  T, 

n*l 


(1.18) 


where  c ,  is  the  congnment  of  s^(t)  along  $  (t), 

N  N 

T 

cN  -  J*  s°(t)tN(t)dt  ;  N  -  1,2,3,...  .  (1.19) 

Confining  equations  (1.18)  and  (1.11),  it  can  be  shown  that  the  equation 
specifying  cj^  is 

T  ® 

J*  |s°(t)|2dt  -  Z  c2(l  +  2ff0AN)(l  +  a0/XN)"2  -  A  .  (1.20) 

0  n-1  ” 

The  least-favorable  signal  of  (1.12)  becomes 

00 

sL(t)  -  Z  c  (l  +  a0/XN)‘SN(t);  0<  t<T  ,  (1.21) 

n-1  w 

and  the  quantity  specifying  the  worst  filter  performance  in  (1.13)  is 
given  by 

T  00 

<hR,sL>  |s°(t)|2dt  -  Z  c2(l  +  a  A  )"1-A]  .  (1.22) 

U  0  n-1  N  U  N 

These  quantities  of  Interest  can  be  determined  when  the  eigenfunctions 
and  eigenvalues  of  the  noise  autocorrelation  are  known;  subsequently,  the 
robust  filter  and  its  performance  will  be  specified  by  equations  (1.18) 
through  (1.22).  The  following  chapters  will  investigate  robust  filter 
performance  when  the  noise  process  tjR^}  Is  determined  to  be:  1)  triangular 
kernel  noise,  2)  ideally  bandlimited  noise,  and  3)  the  Wiener  noise  process. 
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CHAPTER  II 

THE  TRIANGULAR  KERNEL  NOISE  PROCESS 

The  "triangular"  kernel  noise  process  has  been  shown  by  Papoulis  [4]  to 
be  a  model  for  random  binary  transmission.  The  autocorrelation  function  is 
given  by 


.0  if|t-u|>TN 

,l-|t-u|/TN  if|t-u|<TN 


(2.1) 


As  an  example  of  channel  contamination  by  extraneous  signals,  such 

noise  pulses  may  have  pulsewidths  less  than,  equal  to,  or  greater  than  the 

binary  signal  of  interest;  i.e.,  TN  in  (2.1)  is  the  pulsewidth  of  an 

unwanted  random  binary  signal  in  a  communications  channel  and  T  *  [t  i4  T„} 

s  s  N 

represents  the  coincident  interval  of  the  robust  filter  and  the  binary 
signal  of  interest. 

Papoulis  [4]  and  Thomas  [6]  have  calculated  directly  from  the  binary 
pulse  structure  the  autocorrelation  function  of  triangular  noise  processes. 
Their  procedures  are  readily  applied  to  the  case  T  >  TM  for  calculation 
of  {Rjj};  however,  for  a  Mercer  expansion  of  R^(t,u)  on  [O.T^],  the  filter 
solution  of  (1.18)  cannot  be  expressed  in  terms  of  eigenfunctions  and 
eigenvalues  of  £r^}  with  absolute  and  uniform  convergence  assured  on  the 
larger  interval  [0,T  ].  Hence,  the  case  T  >  T„  is  not  considered  in 
this  study. 

Kailath  [5]  examined  in  depth  the  triangular  kernel  for  Tg  <  T^,  as 
an  example  of  kernels  with  nonrational  spectra.  In  particular,  the  power 
spectral  density  function  for  (2.1)  is  given  by 
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S  (f)  -  NqTjj  sinc2(fTN),  (2.2) 

where  Nq  is  the  noise  power  and  f£  -  T^*  is  the  bandwidth.  To  derive  the 
eigenvalues  and  eigenfunctions  of  [r^} ,  define  the  sequence  (PQ>  n  »  1,2,3,...} 
by 


(N-l)rr  <  8  T  <  Ntt  ,  N  -  1,2,3,... 

N  S 


(2.3) 


and 

tan  h  0nT8-±  tTs/(2-Ts)  ][%  P^f1  (2.4) 

where  T  ^  1,  and  thus  0  <  T  <1.  Then  the  eigenvalues  of  (2.1)  are 
N  3 

given  by 


*N 


nq/Pn. 


N 


1,2,3,. 


(2.5) 


and  the  corresponding  orthonormal  eigenfunctions  are 

T 

^(C)  -  <Pn(t)/[J*  |<pn(t)|2dt]*,  0<  ts  Ts, 

N  -  1,2,3, . . .  (2.6) 

where 

c^(t)  -  sin(PNt)  +  cos(PNt);  0  <  t  <  Tg, 

N«  1,2,3,...  (2.7) 

and 

h  -  [T,/(2-T,)H%  V.’'1  • 

Figure  2.1  justifies  the  approximation  implicit  in  the  above  derivation 


P  T  **  $  N/2}  2rr,  N  -  1,2,3,...  (2.8) 

N  9 

where  fx}  denotes  the  greatest  Integer  less  than  or  equal  to  x. 


“■  -  -- 
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Figure  2.1.  Graphical  solution  of  transcendental  equation. 


Straightforward  confutation  using  (2.4)  yields 


I 


dt 


Ta/2[l  +  cos(PNTs)sln(23NT#)/(2aNTt)][l  +  Kjl; 

N  -  1,2,3,... 


(2.9) 


It  then  follows  to  solve  equations  (1.18)  through  (1.22)  for  a  nominal 
signal  s®(t)  and  a  given  degree  of  distortion  A  to  determine  the  associated 
robust  filter  ami  its  performance.  Several  root-finding  techniques  are 
available  to  solve  equation  (2.4);  however,  in  the  case  where  T  /T„  «  1 
(narrow-bandwidth  noise),  the  approximation  in  (2.8)  derived  from  the 
graphical  solution  in  Fig.  2.1  becomes  highly  accurate,  and  the  two 
sequences  and  have  the  following  approximations: 


'  2[Ts(2  -Ts)f%  N-l 

fN/2*  2tt/T  -  (-l)Nt2/(2-T  )  ]  [tN/2*nf 1  ;  N>1 

k  S  8 


(2.10) 


*NoV2-V 


;  N-l 


(2.11) 


[  k  N0[Ts/(fN/2^rr)]  ;  N>1 


As  an  illustration  of  robust  filter  development,  consider  the  nominal 
baseband  signal 


s  (t)  -  EQ  ;  0  <  t  <  T 

E0  »  0  . 


(2. 12) 


From  (1.19)  we  have 


CS  •  Ig  E0  ♦N(t>d,; 

T 

a  E0(l  +  Kj)8inOHT)P;1/[J'  |^(t)2dt]%  N  -  1,2,3 .  (2.13) 

Note  Chat  (2.7)  can  be  rearranged  via  (2.4)  Co  yield 

cpN(t)  *  cos  PN(t -Ta/2)/cos(PNT8/2)  .  (2.14) 

The  robust  filter  Is  then  found  from  (1.18),  specifically 

hR(T-t)  -  hR(t) 

-1  2  " 

-  EQa0A[l  -  (1+KX)  2  yN  cos  PN(t-Ts/2)];  0  <  t  <  Tg  (2.15) 

n«  1 

where 

T 

VN  -  sin(PNT8)(PN(l+a0/XN))"l/J,o|c^(t)|2dt  N-  1,2,3,...  .  (2.16) 

The  robust  filter  performance  for  triangular-kernel  noise  Is  plotted 

In  Figures  2.4  and  2.5.  These  measures  can  be  compared  with  the  filter 

performance  for  the  Gauss-Markov  noise  condition  as  the  distortion  In  the 

signal  structure  Is  Increased.  Results  from  the  Gauss-Markov  study  by 

Poor  [1]  are  plotted  In  Figures  2.2  and  2.3.  An  examination  of  the  quantities 

In  equations  (2.15)  and  (2.16)  reveals  that  the  filter  Impulse  response, 

normalized  to  agE^Si^t),  depends  chiefly  on  parameters  (°q/NqT8)  and 

bandwldth-pulsewldth  product  (TgTN*  s  afTg;  recall  is  defined  as  unity 

In  this  development).  The  distortion  factor  (0q/NqT8)  Is  plotted  versus 

signal  distortion  in  Figure  2.4,  for  the  cases  arT  * 

8 


i 


0.1  and  1.0. 


Mini  max  Error  Probability 


Normalized  Maximum  Distortion,  A/EqT 


F*-«*7* 

Figure  2.3.  Maximin  error  probability  versus  maximum  distortion  for 
narrowband  Gauss -Markov  noise. 


Normalized  Distortion-Noise  Energy  cr0/N0T# 
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Figure  2.4.  Distortion  factor  versus  maximum  signal  distortion  for 
triangular  kernel  noise. 


The  robust  filter  is  calculated  via  (2.15)  for  all  conditions  of 
2  R 

(ofT  )  and  (A/EnT  ).  In  all  cases,  h  (t)  is  symmetric  about  t  ■  T  /2, 

S  U  S  8 

as  predicted  by  (2.14).  The  impulse  response  is  lower  at  this  midpoint 

£ 

than  at  the  endpoints,  but  this  amplitude  differential  in  h  (t)  is  reduced 
2 

as  (A/EqTs)  increases.  This  reduction  of  endpoint  singularities  is 
consistent  with  solution  techniques  for  integral  equations,  where  the 
introduction  of  singularities  into  the  kernel  R^(t,u)  (physically,  a  white 
noise  component  in  the  noise)  reduce  or  eliminate  the  need  for  singularities 

D 

in  a  continuous  square-integrable  solution  such  as  h  (t) . 

To  illustrate  this  concept,  Kailath  [5]  has  shown  that  the  solution  to 
(1.7),  with  (R^(t,u)}  being  the  triangular  kernel,  the  signal  of  Interest 
being  defined  as  in  (2.12),  and  with  A  *  0,  is  an  Impulse  response  consisting 
entirely  of  singularities  at  the  endpoints. 

hR(t)  -  A  6(t)+A  6(t-T  );  OStil  (2.17) 

s  s 

where  A  is  a  positive  constant  determined  by  boundary  conditions.  Thus, 
the  receiver  disregards  the  signal  over  the  entire  observation  Interval 
except  at  the  endpoints  of  the  interval.  This  surprising  result  is  more 
clearly  understood  when  the  triangular  kernel  is  viewed  as  the  model  for 
random  binary  communication.  The  filter  solution  (2.17)  is  then  seen  to 
minimize  the  variance  of  the  noise  process;  i.e.,  this  impulse  response 
achieves  optimal  signal-to-nolse  ratio,  assuming  zero  white  noise 
component  (A  *  0) . 


AC  Che  opposite  extreme,  as  Che  channel  discortlon  Increases  without 

£ 

bound,  Che  slightly  lower  response  of  h  (c)  at  t  *  Tg/2  approaches  the 

maximum  amplitude  at  the  endpoincs.  The  former  singularities  are 

£ 

decreasing  in  amplitude  as  h  (t)  becomes  essentially  a  constant  on  the 

interval  [0,T  ];  i.e.,  h  (t)  approaches  the  matched  filter  solution  for 
s 

detection  in  a  predominantly  white  noise  environment 


hR(t)  -  k  s°(T-t)  «  k  Efl;  0  <  t  <  Tg  , 


(2.18) 


where  k  is  a  scale  factor  irrelevant  to  threshold  detection.  Thus,  for 
both  extremes  of  signal  distortion,  the  robust  filter  computed  by  (2.15) 
is  seen  to  be  the  optimal  impulse  response  for  maximum  signal- co-noise 
ratio  (SNR) . 

The  maximin  error  probability  (Pe>  is  a  filter  performance  measure 

£ 

equivalent  to  worst  case  SNR.  The  robust  filter  max  P  (h  ),  as  given 

s(t)  €«/  6 

by  (1.13)  and  (1.22),  is  plotted  in  Figure  2.5  for  the  case  arT  *  0.1  and 

s 

values  of  SNRq  *  10,  20,  and  30.  SNRq  is  defined  as 


SNR-  -  lim  [<hR,sL)l  -  2  cl/\ 
A  -*  0  n*l 


(2.19) 


Pe(SNRQ)  -  i-0[(SNRor] 


(2.20) 


The  robust  filter  performance  for  the  triangular-kernel  noise, 
measured  by  P ,  is  qualitatively  similar  to  the  results  in  the  Gauss-Markov 
noise  cases.  Performance  is  seriously  degraded  for  small  values  of  channel 
distortion;  this  supports  the  suggestion  by  Poor  [1]  that  this  detection 


procedure  is  very  sensitive  to  I ^  distortion. 
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Computations  to  produce  Figures  2.4  and  2.5  were  performed  by  digital 
computer,  which  required  the  truncation  of  infinite  series  [XN}  and  [$.,(t)}. 

Tw  w 

0  2 

To  minimize  truncation  error,  the  quantity  [J*  |s(t)|  dt]  is  approximated 
N  2  0 

by  E  (c  )  ,  and  (1.22)  becomes 

i  n 

n*l 


<hR,sLV  -  a 


■if j  2 

E  c 

i  n 
U  n**l 


<l+W 


-1 


m 

i 


(2.21) 


For  sufficiently  large  values  of  N,  all  significant  eigenvalues  are  Included 

in  the  computations,  and  the  error  is  inconsequential.  The  estimate  in 

R  L 

(2.21)  is  less  than  the  true  value  of  (h  ,s  )  for  all  values  of  N,  thus 
providing  an  upper  bound  for  probability  of  error.  The  difference  term  for 
the  estimate  is  the  summation 

<hR,sL>-<hR,sL>’  -  E  ^(l  +  X/a  r1  .  (2.22) 

n*N+l  n  U 

This  term  is  made  sufficiently  small  by  choosing  the  proper  value  of  N. 

The  singular  nature  of  the  zero  distortion  filter,  as  seen  in  (2.17), 
causes  the  detection  to  be  sensitive  to  small  amounts  of  distortion. 

Helstrom  [7]  notes  that  this  singularity  can  occur  when  the  noise  auto¬ 
correlation  function  has  no  singular  (or  white  noise)  component.  Such 
sensitivity  should  be  more  pronounced  as  the  noise  processes  under  considera¬ 
tion  become  more  unlike  the  purely  white  noise  case.  A  study  which  assumes 
{R^(t,u)}  to  be  ideally  band  limited  white  noise  provides  a  striking  example 
of  a  robust  filter's  Increasing  sensitivity. 
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If  the  noise  structure  under  consideration  is  to  be  second -order 
stationary,  with  uniform  spectral  density  in  (-a, a)  and  zero  elsewhere,  then 
equation  (3.4)  becomes  the  Karhunen-Loeve  representation  of  that  bandlimited 
white  noise.  This  representation  has  an  expansion  in  terms  of  the 
eigenfunctions  ;  n  -  1,2,3,...}  and  eigenvalues  (X^;  i  •  1,2,3,...} 

as  defined  in  (3.3)  and  (3.4)  on  the  interval  -T/2  <  t  <  T/2. 

The  equation  (3.4)  has  been  Investigated,  and  its  solutions,  called 
prolate  spheroical  wave  functions,  have  been  tabulated.  These  eigenfunctions 
are  seen  to  depend  on  the  product  c  ■  orT/2,  or  2c  *  aT.  Note  that,  since  the 
solutions  are  defined  on  a  fixed  time  interval  [ -T/2, T/2],  a  change  in  the 
constant  c  denotes  a  change  in  the  bandwidth  of  the  noise  process.  It  will 
be  seen  that  a  progressive  study  of  the  aforementioned  singular  nature  of 
a  zero  distortion  filter  is  possible  as  the  robust  filter  performance  is 
determined  for  c  =  0.3,  1.0,  2.0,  and  4.0. 

Numerical  computations  for  this  study,  as  specified  by  (1.18)  through 
(1.22),  will  be  in  terms  of  the  first  four  eigenvalues  and  eigenfunctions 
only,  but  the  properties  of  the  above  series  of  bandlimited  functions  reduce 
the  effects  of  truncation  errors.  In  particular,  we  have  that  the  model 
for  the  signal  class  defined  in  (1.8),  and  the  assumed  noise  process  are 
real  signals,  and  as  such  are  both  timelimited  and  bandlimited.  This 
consideration  is  justified  in  work  [11]  by  Sleplan.  It  follows  from 
Property  1  that  representation  of  real  signals  by  the  bandlimited  functions 
jj(t)}  will  achieve  goodness  of  fit  in  the  interval  [-T/2,T/2];  l.e., 
given  a  member  of  the  signal  class  s(t)  €  J  and 
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with  the  CcN}  given  by  (1.19),  then  the  error  in  the  fit  of  *N(t)  to  s(t)  is 

T/2  • 

-T/2  J*  (s(t)  -sN(t))2dt  -  Z  cjx,,  .  (3.6) 

N+ 1 

The  {X^3  sequence  is  seen  to  approach  zero  rapidly  for  sufficiently  large  N, 
thus  (3.6)  will  be  small  for  such  N.  Accuracy  requirements  and  the  behavior 
of  the  eigenvalues  determine  N.  The  behavior  is  discussed  in  [12],  which 
makes  the  statement  that,  given  a  process  timelimited  to  a  T  second  interval 
and  bandlimited  to  (-a, a),  there  are  only  (oT/tt)  +  1  significant  eigen¬ 
values.  Confirmation  and  more  precise  statements  are  in  [8],  [9],  and  [10]. 
The  highest  value  of  aT  in  this  work  is  eight;  therefore,  the  value  N  *  4 
insures  that  all  significant  eigenvalues  are  included  in  the  computations. 

The  robust  filter  performance  shows  a  departure  from  trends  established 
earlier  by  triangular  and  Gauss-Markov  noise  models  for  the  parameter 
(ct0/NqT).  In  Figure  3.1,  the  normalized  distortion  factor  for  a  given 
fraction  of  signal  distortion  is  seen  to  increase  with  Increasing  product 
(oT).  Note  that,  unlike  the  previous  noise  cases,  the  computations  for  CTq 
in  the  ideal  bandlimited  case  use  a  pulsewidth  which  is  constant  for  all 
values  of  aT.  Thus,  dissimilar  performance  measures  are  anticipated  for 
this  normalized  parameter;  yet  Figure  3.1  is  of  interest  because  the 
equation  (1.20)  defining  Oq  represents  the  imposition  of  an  energy 
restriction  upon  the  signal  of  interest,  resulting  in  a  least  favorable 
signal  as  defined  in  (1.21).  It  is  seen  that,  for  a  given  value  of  arT,  the 
factor  <7q  is  invariant  to  changes  in  Eg;  i.e.,  changes  in  the  energy 


Normalized  Maximum  Distortion,  A/EnT 
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Figure  3.1.  Distortion  factor  versus  maximum  signal  distortion  for 
ideally  bandlimlted  noise. 
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0  2 

|| a  ||  .  It  is  apparent  in  Figure  3.1  that  this  energy  restriction  Cq  varies 

with  the  product  ckT.  More  precisely,  Cq  is  a  function  of  the  eigenvalues 

2 

and  of  the  terms  [c^},  which  represent  the  distribution  of  signal  energy 
in  the  eigenvectors. 

The  increase  in  a  signal  energy  restriction  with  increase  in  arT  should 
affect  maximln  error  probabilities,  and  be  apparent  in  a  comparison  of  P_ 
attained  by  filters  designed  for  noise  models  of  different  c/E.  The 
variation  of  P  vs.  oT  shown  in  Figure  3.2,  however,  reveals  that  any 
effect  on  relative  SNF  by  this  trend  in  Oq/NqT  vs.  oT  is  masked  by  the 
singular  characteristics  of  the  2ero-dis tortion  filter.  Figure  3.2  plots 
error  probabilities  for  SNRq  *  10  and  arT  «  1,  2,  4,  and  8.  The  calculations 
from  SNR  ”  20  and  30  shewed  similar  performance  characteristics,  namely 
that  signal  detection  becomes  increasingly  sensitive  to  small  amounts  of 
channel  distortion  as  arT  decreases  toward  the  narrowband  case.  For  the 
bandlimited  white  noise  process,  a  decrease  in  arT  is  a  measure  of  the 
process'  departure  from  a  purely  singular  autocorrelation;  therefore, 
severe  performance  degradation  is  predicted  as  the  filter  itself  must 
become  increasingly  singular  in  nature. 

B 

The  robust  filter  impulse  response  h  (t)  is  symmetric  about  the  time 
interval  midpoint,  with  evidence  of  suppressed  singularities  on  interval 
endpoints.  There  is  also  evidence  of  local  maxima  and  minima  in  the  inter¬ 
val  of  the  response,  but  the  resolution  of  the  tabulated  data  and  graphs 
of  the  functions  did  not  allow  plotting  of  h  (t)  vs.  time  with 

sufficient  accuracy  for  a  detailed  study  of  impulse  response.  Nor  could 
data  sources  support  closed  form  analytical  expressions  for  the  eigenfunctions 


FP-7244 

Figure  3.2.  MaximLn  error  probability  versus  maximum  distortion  for 
ideally  bandlimited  noise. 
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or  Che  robusc  filter,  prompting  development  of  useful  error  bounds  for 
calculations.  Such  bounds  were  supplied  by:  1)  assurance  that  four 
eigenvalues  provide  sufficient  accuracy  for  the  included  values  of  c*T; 

2)  the  proper  truncation  of  the  eigenvalue  summations  in  the  computations , 
providing  reasonable  lower  and  upper  bounds  for  signal-to-noise  ratio  and 
maximin  error  probability,  respectively.  The  truncation  is  as  described 
in  (2.21)  and  (2.22). 
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CHAPTER  IV 

THE  WIENER  NOISE  PROCESS 

The  Wiener  noise  process  was  developed  as  a  model  for  Brownian  motion. 
It  is  of  interest  because  a  large  class  of  processes  are  generated  from 
the  Wiener  process:  many  others  can  be  transformed  into  a  Wiener  process. 

As  an  example,  a  binary  signal  scheme  is  input  to  an  integrator  before 
transmission;  it  can  be  shown  that  if  the  signal  is  corrupted  by  a  white 
noise  process,  e.g.  thermal  noise,  the  noise  output  of  the  integrator  is 

t 

Y(t)  -  f  N(T)dT  (4.1) 

0 

where  N(t)  is  white  noise  with  covariance 

Rj^Or)  -  ct26(t)  .  (4.2) 

Y(t)  is  a  sample  function  of  a  Wiener-Levy  (or  Wiener)  process  with  the 
properties : 

Y(0)  -  0 

E(Y(t)]  -  0  (4.3) 

E[Y2(t)]  -  o2t  (4.4) 

and  the  density  function 

PYtY<t) ]  -  (2tt  <J2t)"^  exp[-Y2(t)/(2o2t)]  .  (4.5) 

The  Wiener  process,  then,  is  a  nonstationary  random  process.  Using  the 
above  properties,  the  covariance  kernel  is  found  as 
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Ky(t,u) 


2 


a 


min(t,u) 


a2t 


t  ^  u 


2 

\  a  u 


,  U  <  t 


(4.6) 


This  function  will  be  continuous  on  a  given  interval,  thus  the 

2 

Karhunen-Loeve  expansion  exists  for  0  <  t,u  <  T  in  [0,T]  .  The  equation 
of  interest  to  generate  the  series  expansion  (1.15)  becomes 

t  T 

X  *00  «  CT2  J  u^r  (u)du  +  a2t  J  *  (u)du  .  (4.7) 

0  t 

The  solution  sequences  {XN}»  C$N  (c)}  >  end  CcN)  needed  to  represent  the 
robust  filter  for  the  Wiener  noise  process  become 

XN  -  4N0(T/tt)2/(2N-1)2  ;  N  -  1,2,3,...  (4.8) 

i|rN(t)  -  sin(o^t)  ;  0  <  tS  T 

•£  ■  »oAii  •  <4-9> 

The  nominal  signal  is  the  integrator  output 


»0(t)  -  E0t  ;  0  <  t  <  T 


and  the  coefficients  [c^;  N  -  1,2,3,...}  become 

T  -(-1)NEqT3/2  25/2 

*  J  Egt  ^„(t)dt  *  2  2 

N  0  U  N  tt  (2N-ir 


We  also  have  from  (1.18)  that 


OD 

hR(t)  -  Oq1  Eq[  (T-t)  -2T  2  kJJ1  cosOi^t)] 

n»l 


(4.10) 


(4.11) 


(4.12) 
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where 

kjj  -  (I+cJq/^XN -%)2  tt2  .  (4.13) 

O 

The  expressions  for  XN,  (t) ,  Cjj,  and  h  (t)  are  unconplicated,  which 
will  allow  the  robust  filter  performance  to  be  derived  analytically.  An 
inspection  of  (4.12)  and  (4.13)  reveals  that  the  normalized  Impulse  response 

*1  R 

aQ  eq  h  (t)  will  depend  on  the  parameter  (aQ/NQT)  only.  This  distortion 

factor  is  specified  in  (1.20)  and  can  be  expressed  analytically  by  defining 

a  variable  m<  . : 

A 

T 

€[0,0.5] :A  -  ^  J*  | s0(t) |2  dt  ;  (4.14) 

using  the  significant  terms  in  the  sunmation  in  (4.12),  the  following 
approximation  determines  Oq  as  a  function  of  (A/||s^||2) ,  where 

llsV  *  /  |s°(t)  |2  dt  -  Eq  T3/3  ,  (4.15) 

0 

yielding 

a0~  4  N0t2/[tT2(1  -M.^)]  .  (4.16) 

R  L 

The  worst  case  filter  performance  is  specified  by  the  quantity  (h  ,s  ). 

This  expression  of  SNR,  derived  from  (1.21),  becomes 

<hR,sL>  *  oj1  [i  E3  T3(l  -M-a)J  .  (4.17) 

Using  (4.16),  this  performance  measure  becomes 

<hR,sL)  -  tt2  E2  T(l-^)2/(12  Nq)  ; 


(4.18) 
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This  result  can  be  verified  by  considering  the  derivation  for  the  zero 
distortion  case;  i.e.,  using  equation  (2.19),  the  limiting  case  as  (A  *•  0) 

lim  <hR,sL)  «  (c.)2/X.  ;  (4.19) 

A  -  0  L  1 

only  the  first  eigenvalue  is  considered  significant.  This  approximation, 
which  is  derived  independently,  is  practically  identical  to  (4.17)  for 
(^-0): 

<hR,sL>|  -  (tt2/12)  (E2  T/Nq) 

A*0 

-  (c1)2/X1=(8/tt2)(Eq  T/Nq)  .  (4.20) 

The  derivation  (4.17)  has  been  proven  accurate;  however,  the  bounds  for 
the  approximation  error  for  (4.18)  are  a  function  of  the  convergence  of  a 
p-series.  The  series'  first  term  is  proportional  to  the  first  eigenvalue 
and  is  seen  to  produce  error  bounds  too  large  to  allow  confidence  in 
computations  using  (4.17)  or  (4.18);  consequently,  the  analytical  derivation 
is  useful  only  when  general  trends  are  to  be  investigated,  or  in  the  special 
case  where  the  transmitted  signal  s^(t)  is  proportional  to  the  N-th 
eigenfunctions  of  [^(t.u)};  e.g.,  s°(t)  *  a  tN(t).  Poor  [1]  has  shown 
that  the  solution  to  (1.14)  becomes 

hR(t)  -  a  x"1  (I+Oq/Xj,)*1  *N(T-t),  0  <  t  <  T  ;  (4.21) 

which  is  matched  filter  solution  for  coherent  detection,  with  scale  factor 
(I+Oq/X^)"2  .  Such  a  solution  would  have  no  p-serles  involved  in  the 
calculations;  thus,  the  results  would  possess  both  accuracy  and  tight 
error  bounds  based  upon  truncation  of  Insignificant  terms. 
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The  robust  filter  performance  in  noise  modelled  as  a  Wiener  process 
is  demonstrated  in  Figures  4.1  and  4.2.  These  performance  measures  have 
error  bounds  obtained  by  truncation  of  the  insignificant  higher  eigen¬ 
values  (N  >  20).  The  parameter  (Cq/NqT)  is  calculated  as  accurately  as 

possible,  and  plotted  in  Figure  4.1  vs.  [A/||s^||^].  This  term  is  then  used 

R  L 

in  (1.22)  to  calculate  a  lower  bound  for  worst  performance  (h  ,s  ). 

Finally,  this  lower  bound  for  SNR  is  used  in  (1.13)  to  yield  an  upper  bound 
for  the  maximin  error  probability,  which  is  plotted  in  Figure  4.2  for 
values  of  SNRq  *  10,  20,  and  30. 

The  performance  of  this  filter  continues  the  trends  noted  earlier 

for  the  robust  filter.  In  particular,  the  detection  method  appears  sensitive 

to  the  L2  distortion,  as  SNR  is  greatly  reduced  for  small  fractions  of  signal 

R 

structure  distortion.  The  impulse  response  of  h  (t)  is  essentially  that  of 
the  matched  filter  for  a  signal  defined  as  in  (4.10);  however,  the  increase 
in  signal  distortion  produces  a  decrease  in  the  peak  amplitude  of  the 
response.  This  is  similar  to  an  estimation  problem  where  the  robust 
filter  must  operate  on  an  increasingly  distorted  signal  to  obtain  an 
estimate  of  the  signal  waveform  with  which  to  "match",  for  optimum  SNR. 

The  filter's  confidence  in  the  received  waveform  will  diminish  as  distortion 
increases,  and  a  signal  estimate  oust  be  derived  by  the  designed  optimization 
scheme;  e.g.,  minimum  mean-square  error  point.  The  decrease  in  peak 

O 

amplitude  of  a  response  h  (t)  generated  by  such  procedures  thus  yields 
evidence  of  both  the  desired  signal  structure  and  its  increasing  distortion 
by  channel  noise. 


Normalized  Distortion- Noise  Energy  <r0/N0T 


Minimax  Error  Probability 
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Figure  4.2.  MaximLn  error  probability  versus  maximum  signal  distortion 
for  the  Wiener  noise  process . 


CHAPTER  V 


SUMMARY  AMD  SUGGESTIONS 


Matched  filtering  is  seen  to  be  sensitive  to  signal  distortion  when 

designed  for  three  noise  structures:  the  triangular-kernel  noise,  ideally 

bandlimited  noise,  and  Wiener  noise  processes.  It  should  be  noted  that 

R  L 

pessimistic  scenarios  have  been  presented:  1)  the  quantity  (h  ,s  )  specifying 
the  worst  performance  was  used  to  calculate  signal-to-noise  ratios;  2)  in  all 
cases,  an  upper  bound  for  maximin  error  probability  was  plotted  as  the 
actual  filter  performance. 

If  the  robust  filter  detection  procedure  is  to  be  pursued,  it  would  be 
advantageous  to  consider  l^-norm  distortion  signal  models  such  as  (1.8) 
which  limit  the  effective  maximum  A*  Successful  research  in  radar  pulse 
design  is  an  example  of  improved  signal  structure  modelling. 

In  addition,  future  noise  models  to  be  considered  might  feature  a 
singular  (white  noise)  component  in  the  noise  autocorrelation  function; 
e.g.,  if  the  noise  structure  is  ideally  bandlimited  as  in  (3.1),  define  a 
new  noise  process 


Sy(w)  -  Sx(w)+N0/2 


'  V2+s0’  M  <0f 

,  N0/2  ,  |w|  >  a  . 


(5.1) 


This  technique  may  reduce  the  robust  filter  sensitivity  to  distortion, 
although  the  Mercer  expansion  of  R^.(t,u)  or  the  solution  to  the 
Karhunen-Loeve  expansion  might  become  difficult  to  determine. 
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